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Abstract. Recently, Holz & Wald have presented a new 
method for determining gravitational lensing effects on, 
e.g., supernova luminosity versus redshift measurements 
in inhomogeneous universes. In this paper, their method is 
generalized in several ways: First, the matter content is al- 
lowed to consist of several different types of fluids, possibly 
with non- vanishing pressure. Second, besides lensing by 
simple point masses and singular isothermal spheres, the 
more realistic halo dark matter distribution proposed by 
Navarro, Frenk & White (NFW), based on N-body simu- 
lation results, is treated. We discuss various aspects of the 
accuracy of the method, such as luminosity corrections, 
and statistics, for multiple images. We find in agreement 
with other recent work that a large sample of supernovae 
at large redshift could be used to extract gross features 
of the mass distribution of the lensing dark matter halos, 
such as the existence of a large number of point-like ob- 
jects. The results for the isothermal sphere and the NFW 
model are, however, very similar if normalized to the ob- 
served luminosity distribution of galaxies. We give conve- 
nient analytical fitting formulas for our computed lensing 
probabilites as a function of magnification, for several red- 
shifts. 

Key words: cosmology: theory, gravitational lensing, 
dark matter - galaxies: halos 



1. Introduction 

Ever since the foundation of modern cosmology based on 
general relativity, there have been attempts to determine 
the cosmological parameters of our universe. Gravitational 
lensing has been recognized as one of several such tools. 
For example, the results of radio surveys for gravitational 
lenses have recently been used to place significant con- 
straints on cosmological models (Falco et al. 1998 ). It has 
also been realised that observations of supernovae at high 
redshift can be used for this purpose (Goobar & Perlmut- 
ter 



cosmological constant. In fact, several collaborations with 
this in mind arc in progress, and the first sets of data 
show an intriguing hint of a non-vanishing cosmological 
constant in a universe consistent with having a flat geom- 
etry (see, e.g., Riess et al. 1998; Perlmutter et al. 1999). 



Although the two groups which have published results mu- 
tually agree on the best-fit parameters, it is important to 
note that the effects of geometry are small (on the order 
of half a magnitude), and the need to go to even higher 
redshift to get larger effects is obvious. When observing 
such distant sources, at redshift greater than unity, it is 
necessary to estimate the effects of lensing due to inhomo- 
geneities in the matter distribution. This will be of grow- 
ing importance when supernovae at very high redshifts 
become accessible, e.g., with NGST, the Next G eneration 
Space Tel escop e (Miralda-Escude & Rees 1997 ; Dahlen & 
Fransson 1999] ) or with the dedicated supernova search 
satellite, SNAP (Perlmutter et al., private communica- 
tion). 

The literature on gravitational lensing is quite rich. 
Much of the history and developments up to the early 
1990s can be found in the excellent textbook by Schnei- 



der, Ehlers & Falco (1992). A new method for examin- 
ing lensing effects has recently been proposed by Holz & 
Wald (1998; HW), one which has the virtue of lending it- 
self easily to numerical calculations. In that work it is also 
shown that given a large and deep enough sample of stan- 
dard candle supernovae, one may in principle distinguish 
between different matter distributions in the galactic ha- 
los responsible for the lensing. In particular the case of 
point masses gives a quite different distribution of mag- 
nifications (or de-magnifications) compared to the case of 



singular isothermal spheres (see also Metcalf & Silk 1999) 



1995), in particular for determining the value of the 
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The method of HW can be summarized as follows: 
First, a Friedmann-Lcmaitre (FL) background geometry 
is selected. Inhomogeneities are accounted for by speci- 
fying a matter distribution in cells which have an aver- 
age energy density equal to that of the underlying FL 
model. Following the history of observed light from a dis- 
tant source, a light ray is traced backwards to the desired 
redshift by being sent through a series of cells, each time 
with a randomly selected impact parameter with respect 
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to the matter distribution in the cell. Between each cell, 
the FL background is used to update the scale factor and 
the expansion rate. By using Monte Carlo techniques to 
trace a large number of such light rays, statistics for the 
apparent luminosity of an ensemble of sources at a given 
redshift is obtained. 

The purpose of this paper is to generalize the method 
of HW in a number of ways and to investigate some issues 
related to the different lensing signatures in the luminos- 
ity distribution of supernovae caused by different mass 
distributions of the intervening lens population. The out- 
line is as follows. In Sect. || it is demonstrated how the 
HW method can be generalized from matter in dust form 
to perfect fluids with pressure. Sect. || is concerned with 
the interpretation of results as probability distributions. 
Sect. |H presents the geodesic deviation and in Sect. ||, we 
summarize and discuss some of the conceptual issues in the 
analysis of HW. In Sect. ||, we include the possibility of 
using the more realistic halo matter distribution proposed 



by Navarro et al. (1997; NFW) based on their N-body sim- 
ulation results, and in Sect. ^ we discuss the mass distribu- 
tion of lensing galaxies. By using realistic density profiles 
and mass distributions of dark matter halos, we can ob- 
tain high accuracy results on gravitational lensing without 
the need of using the extensive full data sets of N-body 
simu lations, as in so me oth er method s (Wa mbsganss et al. 
199S| ; Premadi et al. |l998[ Jain et al. |l999|) . One strength 



of our ray-tracing Monte Carlo method is further that 
it can be continuously refined as more observational in- 
formation, e.g., on galaxy distributions is obtained. Also, 
effects such as absorption by dust and other possibly z- 
dependent effects can be straightforwardly added to the 
algorithm. In Sect. |[ some results from the generaliza- 
tions of the method of HW are presented as luminosity 
distributions. In Sect. ^, we investigate some limitations 
of the numerical implementation of the method and in 
Sect, [l^, we discuss multiple images. Sect. |ll| investigates 
different consistency checks, Sect. |l2] contains convenient 
analytical fits to our numerical results and the paper is 
concluded with a discussion in Sect. O. 



2. Cosmological model 

The starting point of HW is a Newtonianly perturbed FL 
universe. The line element of Robertson- Walker form is 
then 



ds 2 = -(1 + 2<p)dt 2 + (1 - 2cf))a 2 (t)dl 2 k , 



(1) 



where a is the scale factor, and dl\ is a constant-curvature 
three-geometry, for which k 6 {0, ±1} indicates the sign 
of the curvature. (We use geometrized units such that 
G = c = 1). The perturbation <f) is assumed to satisfy 
a number of properties that will be employed when deter- 
mining dominant terms. Firstly, \<p\ <C 1, reflecting that 
<fi indeed is a perturbation. Secondly, inhomogeneities are 



assumed to be significant, so that d(f>/dt is small com- 
pared to spatial derivatives. Thirdly, second-order spatial 
derivatives will be assumed to dominate over products of 
first-order ones (e.g., (d x cj)) 2 <§C \d xx (j)\). In this paper, we 
have generalized the treatment of HW so that the mat- 
ter content of the model can be taken to be a number of 
perfect fluids: 



Tab = PiU a Ub + Pi(u a Ub + 9ab), 



(2) 



where T a ^ is the stress-energy tensor, pi is the energy den- 
sity, pi the pressure of fluid i, u a is the (normalized) fluid 
four-velocity, and g a b is the metric corresponding to the 
line element (|l|). Cosmologically reasonable equations of 
state are linear barotropic ones: pi = oupi. In this paper, 
the cti's are constants with —1 < a. L < 0. However, it 
is straightforward to allow for the equations of state to 
evolve with redshift, i.e., ai = ai(z). 

HW considered dust models with a cosmological con- 
stant, which corresponds to a two-fluid model, the first 
fluid having ai = 0, and the second fluid homogeneously 
distributed with ct2 = — 1 and A = 8np2. The field equa- 
tions, corresponding to Eqs. (HW7,8), become 
The Friedmann equation: 



3H 2 = 8irJ2i 



3k 



2V 



The Raychaudhuri equation: 



-3qH 2 = -Aur^iPi + 3 Pi ) + V 2 0, 



(3) 



(4) 



where the definitions of the Hubble scalar H and decel- 
eration parameter^ q have been made explicit, and where 
V 2 denotes the spatial projection of contracted covariant 
derivatives. The pressure does not enter the Friedmann 
equation, i.e., the Hubble scalar is entirely given by the en- 
ergy densities and the curvature. The deceleration, on the 
other hand, is affected by the presence of pressure. Com- 
paring the above field equations with the field equations 
for the corresponding FL background (with averaged den- 
sities and pressures pi, pi) it can be shown [Eq. (HW11)] 
that <j) satisfies a Poisson equation V 2 (f> — 4:irJ2i°~Pi, 
where Spi = pi — pi, while dpi = Pi — pi vanishes. Thus, 
the pressures result in homogeneous contributions, and 
the specification of an equation of state only makes sense 
in terms of the averaged energy density: Pi = pi = QiPi. 

From the field equations for the background model fol- 
lows 

The energy conservation equation: 



\p i + 3{p i +p i )H = Q. 



(5) 



The name is of historical origin. Actually, according to the 
current supernova cosmology data the expansion of the uni- 
verse is presently accelerating, i.e. q < 0. 
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If the fluids are assumed to couple only to gravity, they will 
be separately conserved, and each component in the sum 
of Eq. (||) vanishes. Note that it is only for dust (p = 0) 
that the mass within a volume co-moving with the Hubble 
flow is constant (p oc a~ 3 ). For other types of matter, the 
mass varies with time because of the work done by the 
pressure. Also note that p = —p results in a constant 
energy density. Thus, a homogeneously distributed fluid 
with p = —p is equivalent to a cosmological constant, or 
in terms of quantum field models, to vacuum energy. It is 
convenient to introduce density parameters 



(6) 



which depend on time and therefore on redshift. Assum- 
ing the above discussed equations of state the time de- 
pendence is governed by the following set of conservation 



equations (see, e.g., Goliath & Ellis 1999) 
Cli = [2q - (3ai + 1)] VtiH. 



HW proceed by performing a number of coordinate 
transformations to co-moving coordinates on isotropic 
form, Eq. (HW16). The inclusion of pressure results in 
a modification of this line element: 



ds 2 = - 



1 + 2$ + AttR 2 



dT 2 



+ [1 - 2$] (dX 2 + dY 2 + dZ 2 ), (8) 
where R 2 w X 2 + Y 2 + Z 2 , and the effective potential 
2ttR 2 



$ = 



satisfies the Poisson equation 

V 2 $ = 471-^/),;. 



(9) 



(10) 



The lookback time of an event at redshift z\ is the time 
difference (with respect to the coordinate t) between the 
event and the present. By using 1 + z = ao/a, we get a 
relation between the lookback time and the redshift for a 
FL universe: 

rao da f Zl dz 



At 



(11) 



To obtain H(z), the Friedmann equation, Eq. (g), is ex- 
pressed in terms of redshift and the present densities poi- 
To do so, first note that the energy conservation equation, 
Eq. (||) , together with the equation of state gives 



| = p oi (l + z)^ 1+ ^. 

The Hubble scalar can then be expressed 



-3(1 + Qi) 



H(z) = H Y, + ^) 3(1+ai) + *M1 + z) 2 



(12) 



(13) 



where we introduced the curvature density parameter 

k 



Ok 



H 2 n 2 ' 
T1 o a o 



(14) 



The luminosity distance to the event at redshift Z\ is given 
by a similar expression (see, e.g., Bergstrom & Goobar 
19991) 



d-L = a (l + z) 



/f 1 f 

J \aoJo 



dz 



sin(a;), k = 1 
fix) = { x, k = 0, 

sinh(a;), k = — 1 



(15) 
(16) 



Note that, when k = ±1, the present scale factor can be 
expressed [using (fil)] as a = (H Vl^ofcl) -1 - 



(7) 3. Area vs. magnification and statistical weighting 



Consider a narrow (strictly speaking; infinitesimal) ray 
bundle connecting a source and an observer as depicted in 
Fig. |l[ The ray bundle can be focused either at the source 
or the observer. The relevant quantity when determining 
luminosity distances is the ratio of the area of a beam 
and the solid angle it subtends. Of course, the smaller the 
distance, the brighter the source. 




Source (observer) 



Observers aperture (extended source) 




Fig. 1. Geometry of ray bundle. 



In the method of HW, a beam focused at the observer 
with infinitesimal <5f2 is integrated backwards using the 
geodesic deviation equation, and the resulting area at the 
source sphere is considered. In the lower panel of Fig. [j], 
the beam has undergone additional focusing. A finite aper- 
ture at the observer will detect light from a larger S£l of an 
isotropically emitting source, and a finite source will sub- 
tend a larger solid angle on the sky. Since surface bright- 
ness is conserved, the source will be magnified. The mag- 
nification, /i, is defined as the ratio of the observed flux 
and the flux of the corresponding image in absence of any 
additional focusing. 



(17) 
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where <5f2o denotes the solid angle subtended by the source 
in absence of lenses. Hence the magnification will be pro- 
portional to the inverse of the area of the beam. In absence 
of any lenses, the ray bundle will have an area correspond- 
ing to the so called empty-beam area (see Sect. [|), and we 
will therefore define the magnification as 



.4, 



empty 

~A~ ' 



(18) 



For further reference, see Schneider, Ehlers & Falco 

As noted in HW, some care is needed when interpreting 
results as probability distributions. This is due to the fact 
that individual ray bundles do not correspond to random 
source positions. This is illustrated in the left panel of 
Fig. U where we consider the situation of a ray bundle 
originating from an observer at z = at the center of a 
source sphere at some fixed redshift. 





Observers sphere 
Fig. 2. Different ray bundles with constant <5£1 



From the right panel of Fig. g we see that focused 
beams sample a smaller fraction of the source sphere. 

Since the fraction of the source sphere being sam- 
pled by a beam is proportional to the area of the beam, 
i.e., inversely proportional to the magnification, magnified 
sources will be overrepresented if not compensated for. In 
HW, this is done by convoluting the obtained probabil- 
ity distribution, no(A), with a factor proportional to A, 
thus obtaining a new distribution n{A) oc A rig (A). The 
method we use is very similar in that we designate a prob- 
ability, proportional to the area of the ray bundle, to keep 
every event. The advantage of this rejection method is 
that each selected event represents a statistically "real" 
event. 

Further discussions of the statistical weight given to 
random line-of-sights can be found in Ehlers & Schneider 

HP - 

4. Geodesic deviation 

The g eodes ic deviation equation for a light ray is (see, e.g. 
Wald 

llX 2 



1984) 



(20) 



where £ a is the deviation vector and k a is a null vec- 
tor tangent to the light ray. The affine parameter A is 
a parametrization along the light ray given by uidX = dt, 
where u> is the frequency of the photons in the beam. It is 
instructive to calculate the deviation in an exact (homo- 
geneous) FL universe, analogous to Eq. (HW21): 



d 2 £ a 
U 2 



= -4™ 2 5> 4 + Pl )r- 



(21) 



Thus, the pressure will cause additional focusing of the 
light beam. Note that a cosmological constant, which cor- 
responds to p = —p, does not affect the deviation. 

To derive the deviation experienced when propagating 
through a cell in the inhomogeneous model, HW start by 
writing the deviation equation in terms of a 2 x 2 matrix 
A a p, such that 



d 2 A a 



d\ 2 



(22) 



where a,0 G {X, Y} (HW28). The properties of the 
geodesic deviation are contained in -4 a /3.[] For example, 
the area (with orientation) of the beam is given by det A. 
In the case of the homogeneous background Friedmann- 
Lemaitre cosmology, we denote the corresponding refer- 
ence area Afl, which will of course depend on the red- 
shift, but also on the cosmological parameters of the FL 
background model. Also, the area of a light beam that 
propagates without encountering any matter will be de- 
noted A emp t y . Sometimes Afl and A em pty are referred to 
as the filled-beam and empty-beam areas, respectively. 

The deviation equation, Eq. (p2|), can now be written 
as a set of difference equations (HW35,36) 



{A a p)x = (A a p) + AX(dA a /dX) o , 

{dA a p /dX) 1 = {dA a (3/dX)o-uJ%(A^) , 

with 

1 



J" 







dXA' 



(23) 
(24) 

(25) 



where the integration is along a straight line through the 
cell. Next, A a p is calculated from the line element Eq. (pj) 



ia/3 : 



-Lo 2 {2d a dp<$> 



Z a d a (Z b d b ^>) + A7rJ2p t 



?M>,(26) 



where rj a b is the Minkowski metric associated with the co- 
ordinates (T, X, Y, Z). Thus, with non- vanishing pressure, 
there will be an additional contribution in J"g that takes 
the form 



rot 




(27) 



2 This matrix is conventionally discussed in con nectio n with 
conjugate p oints for null geodesies, (see, e.g., Wald ', 



(19) ing fc Ellis |l973l) 



1984 



Hawk- 
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with 



J p . — Airupi / dX 



4npi / dZ 




where b is an impact parameter, and R c and V c are the ra- 
dius and volume of the cell, respectively. Assuming linear 
barotropic equations of state, pi = ctipi, Eq. (p8|) becomes 



d O 




(29) 



where Mi is the mass of matter component i in the cell. 

Let us now examine the above result for the filled- 
beam case. Eq. (HW39) gives the quantity J a p for matter 
uniformly distributed in a ball of radius d at the center 
of a cell of radius R c . Setting d = R c in this equation, 
and adding our derived pressure contribution ^2 i J Pi S a /3 
to this yields 



(j fl )*x = (j fl )^ = E( 1 +^)^Wi 



b 2 

m 1 



(J^fy = (J^) Y X = 0, 



(30) 



from which it is seen that a cosmological constant (a 
1 ) does not contribute, in consistency with Eq. (pi|). 
As an interesting example, consider a flat cosmological 
model in which ordinary matter (a = 0) constitutes only 
10 % of the total matter content, i.e., VLm = 0.1, while 
the rest is made up of some homogeneously distributed 
"X-matter" with fix = 0.9. The ordinary matter is as- 
sumed to be in the form of point masses (we will return 
later to more realistic mass distributions). For this cos- 
mology, Garnavich et al. ( 1998 ) have put an upper limit 
on the parameter of the equation of state for an X-matter 
component using Type la supernova data. Their limit is 
ax < —0.60 at 95 % confidence level. Similarl y, the Su- 
pernova Cosmology Project (Perlmutter et al. 1999) has 



deduced a limit of ax < —0.44 (95 % cl) for any value of 
Q.M- 

In Fig. ||, we compare the cases of ax = —1, corre- 
sponding to a cosmological constant, ax = —2/3, which 
may either arise from "quintessence" (i.e., a slowly evolv- 
ing scalar field) (Wang et al. 1999| ) or from a domain wall 
network (Battye et al. 1999| ). For reference, we also show 
the result for ax — 0, corresponding to dust (i.e., equiv- 
alent to an Einstein-de Sitter model with Qm = 1, of 
which 90 % is homogeneously distributed). The choice of 
JIM =0.1 is made to allow for a direct comparison of the 



2 

A 

Afl 
1.5 



(28) o.5 
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Cumulative distribution 

Fig. 3. Beam areas normalized to the filled-beam value of 
the underlying Friedmann-Lemaitre model ^4fl for 10 000 
sources at redshift z = 3 in a Qm = 0.1, Qx = 0.9 uni- 
verse. The matter is assumed to be in the form of point 
masses, and results for the equation of state coefficient 
ax of the homogeneous X component having the value 
(solid line), —2/3 (long-dashed line) and —1 (short-dashed 
line) are shown. The horizontal scale gives cumulative per- 
centage of beams having an area less than a given value. 
Negative areas correspond to secondary images. 

curve corresponding to ax = — 1 in Fig. |^ with Fig. 7 in 
HW, with which it agrees well. 

Note that the normalization differs between the three 
cases since Afl depends on the cosmology. Nevertheless, 
from the slope of the curves, one can conclude that lens- 
ing effects become more prominent for large negative ax- 
This is not what one would expect from studying the 
geodesic deviation equation, Eq. (pil|), since negative pres- 
sure causes less focusing of a beam. However, for negative 
ax, distances will get larger for a fixed redshift, which 
means a larger number of lenses between the source and 
the observer (see, e.g., Turner 1990| ). 



5. Summary of the properties of the lensing 
model 



The method of Holz & Wald ( 1998Q which we adopt has a 
number of properties, which we summarize and comment 
upon here: 

Local nature of lensing 

If 4> is small within distances much smaller than the 
Hubble radius Rh = H , then the curvature (and 
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thereby the lensing effects) is determined locally within 
Rh- Likewise, if we assume that there is no strong 
correlation on scales greater than a clustering scale 
R, it follows that the curvature is determined locally 
within R. 
— Relevant scales 

The question now is what relevant clustering scale to 
use. HW find it convenient to employ the typical sep- 
aration between galaxies, and argue that galaxy clus- 
tering should have a negligible effect on magnification, 
as shown in some of their simulations. They also argue 
that the smallest relevant scale for matter clumping is 
of the order 1O _3 M0, since objects smaller than that 
would not give a coherent lensing effect over the angu- 
lar size of most of the interesting sources. 
Generality 

For the case of point masses, HW demonstrate that 
their results are fairly insensitive to the size of individ- 
ual point masses. Furthermore, their results are rather 
insensitive to clustering. This is tested by cutting out a 
tube in a uniform matter distribution and distributing 
"stars" within this tube. 
Largest effects 

Based upon their results, HW conjecture that random- 
ly distributed point masses should cause the largest 
lensing effects of all possible mass distributions for 
a given average mass density. Note that the point- 
mass case corresponds to the limit in which all mat- 
ter clumps lie well within their Einstein radii [see 
Eq. ©]. 

Comment on relevant scales: 

The Einstein radius is a characteristic length scale in the 
lens plane given by 



Co 



2i?s 



D d D. 



chw " 



d^ds 

D s 



(31) 



where i?schw = 2M is the Schwarzschild radius of the 
lens and Dd s , Dd and D s are angular diameter distances 
between lens and source, observer and lens, and observer 
and source, respectively. The corresponding length scale 
in the source plane, 770, is 



D s DsD ds 
-Jo = ^2R Schw ^- 



(32) 



For lensing to be important, the source size, I, has to 
be smaller than the length scale 770 introduced by a point 
mass lens, I < t/q. The source size of a type I supernova 
is ~ 10 15 cm at maximum luminosity. Putting the source 
at z = 1 and the lens at z = 0.5 in a flat universe with 
&m = 0.3 and the dimcnsionless Hubble parameter h = 
0.65, we obtain 



M > 10" 4 M 



0) 



(33) 



i.e., only compact objects with M > 10 4 M Q should be 
of importance through microlcnsing of supernovae. 



6. The Navarro- Frenk- White distribution 

The work of HW has as a first treatment been concerned 
with the J a p expressions for idealised cases like point 
sources; singular, truncated isothermal spheres (SIS); uni- 
form spheres; and uniform cylinders. However, another 
often-used matter distribution is the one based on the 
results of detailed N-body sim ulatio ns of structure for- 
mation given by Navarro et al. ( 1997 ). They parametrize 
the density profile by 



p(r) 



PcritOc 



(r/R s p[l + (r/R s )VS](P-^S' 



(34) 



where p C rit = 3H 2 /8ir is the critical density, 8 C is a di- 
mcnsionless density parameter and R s is a characteristic 
radius. 

NFW found that halos found in N-body simulations 
ranging in mass from dwarf galaxies to rich galaxy clusters 
can be fit by Eq. (|J) with 7 = 1, /? = 3, 5 = 1. The 
potential for this density profile is given by 



$(r) = -47r/3 crit ^c-Rs — + ^ + const., 



(35) 



where x — r/R s . The matrix J a p can then be obtained 
analytically, see Appendix A. 

The mass inside radius r of a NFW halo is given by 



M(r) = 4Trp crit 6 c R 3 s 



ln(l + x) 



l + x 



(36) 



Since we want to have the average mass density in each 
cell equal to that of the background cosmology, we get the 
relation 



4"7T 

m = —n MPcrit R 3 c . 

Setting M — M(R C ), we obtain 







M 



ln(l + x c ) - 



l+x c 



(37) 



(38) 



where x c = R c /R s . That is, for a given mass M, S c is a 
function of R s . From the numerical simulations of NFW 
we also get a relation between 8 C and R s . This relation is 
computed numerically by a slight modification of a For- 
tran routine kindly supplied by Julio Navarro. Of course, 
one wants to find an R s which gives both the required 
average density in each cell and compatibility with the 
numerical results of NFW. 

Generally, R s will be a function of mass M, the Hubble 
parameter h, Qm, and z. However, we use the result 



from Del Popolo (|1999)) and Bullock et al. (|1999j) that R s 
is approximately constant with redshift. 

A comparison of luminosity distributions between 
point masses, isothermal spheres (SIS) and the NFW mat- 
ter distribution is presented in Sect. 



Lars Bergstrom et al.: Lensing effects in an inhomogeneous universe 



7 



7. Mass distribution 



As discussed above, HW have shown that when studying 
luminosity distributions and other quantities where the 
null geodesies only differ by infinitesimal amounts, their 
method is insensitive to the individual masses and cluster- 
ing of point mass lenses. For point mass lenses, the mass 
distribution and number density will therefore be of small 
importance. However, realistic modeling of galaxies does 
call for realistic mass distributions and number densities, 
i.e., one has to allow for the possibility of R c to reflect the 
actual distances between galaxies, and for galaxy masses 
to vary in agreement with observations. 

One of the advantages of the method of HW is that any 
mass distribution and number density, including possible 
redshift dependencies, can easily be implemented and used 
as long as the average energy density in each cell agrees 
with the underlying FL model. 

To exemplify this, we derive a galaxy mass distribu- 
tion, dn/dM , by combi ning t he Schechter luminosity func- 



tion (see, e.g., Peebles |1993| , Eq. 5.129) 



dn 
0(2/) 



(f>{y)dy, 



(39) 



y = t~ 



L 



with the mass-to- luminosity ratio (see, e.g., Peebles 1993 , 
Eq. 3.39) 

— cx Ml (40) 
L 

Normalizing to a "characteristic" galaxy with L = L* and 



y 



M 
M 

W* 

Using Eq 
dn 

dM K 
S = 



M*, this can be written 

1/(1-/3) 



V s e~ 



j), we find that 



1-P 



(41) 

(42) 
(43) 



In Fig. ^, mass distributions for various values of a are 
displayed. Assuming that the entire mass of the universe 
resides in galaxy halosf] we can write 

J/max 

n(y)M(y)dy = p m . (44) 



In our model with cells joining each other along the paths 
of rays, this means that the underdensity in between galaxies 
is parametrized by taking the halo mass distribution to be 
valid out to large radii. A discontinuity of mass density then 
occurs at the boundary between cells. This has a negligible 
effect on the lensing properties, however, since the mass density 
at the discontinuity is very small. A similar, and in fact larger 
discontinuity occurs for the truncated isothermal sphere where 
the density suddenly drops from typically p ~ 10p cr it to p = 
at the truncation radius. 
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Fig. 4. dn/dM in arbitrary units for f3 = 0.2 in the in- 
terval 0.045 < M/M* < 22.4. The full curve corresponds 
to a = —0.70; the long-dashed to a = —1.07; and the 
short-dashed to a = —1.25. 



Using the Schechter luminosity function and the mass-to- 
luminosity ratio we get 



AL 



^MPcrit 



vdy 



(45) 



Thus, by supplying values for n*, reasonably well-det- 
ermined by observations, and y mm and y max , on which 
the dependence of M* is weak, together with parameters 
a and /3 we can obtain an M* consistent with Qm- 

Note that assuming that the Schechter luminosity 
function and the mass-to-luminosity ratio is valid over 
the entire luminosity range (i.e., putting y m i n = and 
2/max = oo), we get 



AL 



^A/Pcrit 
n*r(p+ 1)' 



(46) 



where p = a + 1/(1 — /?), p > —1. 

A compilation of parameter values can be found in 
[Appendix B[ For 0.5 < y < 2, a = -0.7 and (3 = 0.2, 
Eq. ([46|) can be written 



AL 



7.5- io 13 n 



0.65 



1.9 • 10- 2 h 3 
Mpc 3 



Mr. 



(47) 



7.1. Truncation radii for SIS-lenses 



The SIS model is built on the assumption that stars (and 
other mass components) behave like particles in an ideal 
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gas, confined by their spherically symmetric gravitational 
potential. The density profile is given by 



2tt r 2 



(48) 



where v is the line-of-sight velocity dispersion of the mass 
particles. [Note that this corresponds to Eq. (|34|) with (3 = 
7 = 2, S ?t and v 2 /(2ir) = p alt 5 c R 2 .] The mass of a SIS 
halo truncated at radius d is then given by 



M{d) = / p (r)dV = 2v z d. 
Jo 



(49) 



We want this to be equal to the mass given by the 
Schcchter distribution M, 



2v 2 d = M 



d = 



M* 
2v 2 



M 



(50) 



where, in addition to M* , we have introduced a character- 
istic velocity dispersion v* . Combining the Faber- Jackson 
relation 



— = V 



(51) 



where y is defined in Eq. (|3^) , with the mass-to- luminosity 
ratio, Eq. (flo|), we can substitute for v in Eq. (^o|), 



Mi 

2v 2 



M 



l-2A(l-/3) 



(52) 



For values of i/» and A, see Appendix B . Using Eq. (tf7|), 
we can write the truncation radius for a halo with mass 
M = AL as 



d~ 3.3 ft 



M 



0.65 



1.9- 10- 2 h 3 
Mpc 3 



Mpc. 



(53) 



A comparison between point masses and isothermal 
spheres for sources at z = 3 in a CIm = 1, fl\ = universe 
is presented in Fig. ||, where the full line corresponds to 
matter in the form of point masses. This plot agrees well 
with Fig. 5 of HW. The long-dashed line corresponds to 
matter clumped in isothermal spheres with cut-off radius 
d = 200 kpc. We choose this fix d in order to be able 
to compare directly with Fig. 11 of HW. Both lines agree 
very well with the results in HW. Finally, the short-dashed 
line corresponds to isothermal spheres for which the cut- 
off radius d has been calculated according to Eq. (|52| ) (see 
Sect. 7_A). For the presented model, this results in a typical 
cut-off of d ~ 3 Mpc, which explains why this case comes 
closer to the filled-beam value. 

8. Results 



In Figs. | [14, we compare the luminosity distributions ob- 
tained for NFW halos with the point-mass and SIS mat- 
ter distributions. Comparing Fig. ^| with Fig. 18 and 20 
of HW, we see that the point-mass cases agree very well. 
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Fig. 5. Beam areas normalized to Afl for 10 000 sources 
at redshift z = 3 in a f2j\/ = 1, f^A = universe. Negative 
areas correspond to secondary images. The dashed line at 
the top indicates the empty-beam value A cmpt y. 

The SIS luminosity distributions differs in that our re- 
sults is shifted towards the filled beam value. This is due 
to the fact that the truncation radii used [see Eq. (|^)] are 
an order of magnitude larger than the (constant) value of 
d = 200 kpc used in HW. 

Lensing effects are governed by the lensing potential of 
the galaxy model used, and the number of lenses between 
the source and the observer. As noted in the discussion 
of Fig. [|, a larger number of lenses will tend to make 
lensing effects for a given galaxy model more prominent. 
For a fixed source redshift and a fixed number density 
of lenses, the number of lenses will only be a function of 
the cosmology. We get the largest number of lenses in the 
Om = 0.3, = 0.7 cosmology, followed by the Qm = 
0.3, = case, and the smallest number of lenses for 
Q M = 1, A = 0. 

For point masses, the lensing potential will be unaf- 
fected by the cosmology. As discussed in Sect. ||, point 
masses cause the largest lensing effects for a given average 
mass density. 

For a SIS halo, the truncation radius, d, is proportional 
to VLm [see Eq. ([45]) and Eq. (|52|)]. Since the potential for 
a truncated SIS halo is equal to the point-mass potential 
for r > d (see Fig. |l^), the SIS luminosity distribution 
will be shifted towards the point-mass distribution when 
lowering fijvf. In fact, the agreement between our results 
and those in HW is better for lower JIm, as can be seen by 
comparing the curve corresponding to the SIS luminosity 
distribution in Fig. with Fig. 22 of HW. Taking also the 
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rrac(filled beam) — mag(lensed) mag(filled beam) — mag(lensed) 



Fig. 6. Luminosity distributions for 10 000 sources at red- 
shift z = 1 in a Q,m = 1, = universe. The magni- 
fication zero point is the fillcd-beam value. The full line 
corresponds to the point-mass case; the dashed line is the 
distribution for SIS halos, and the dotted line is the NFW 
case. This plot can be compared with Figs. 18 and 20 of 
HW. 



Fig. 7. Luminosity distributions for 10 000 sources at red- 
shift z = 1 in a CIm — 0.3, J7a = 0.7 universe. The mag- 
nification zero point is the filled-beam value. The full line 
corresponds to the point-mass case; the dashed line is the 
distribution for SIS halos, and the dotted line is the NFW 
case. This plot can be compared with Fig. 22 of HW. 



number of lenses into consideration, we expect the largest 
lensing effects for the J7m = 0.3, Cl\ = 0.7 cosmology, 
followed by the Hm = 0.3, 57a = case, and the smallest 
effects in the f^,/ = 1, = cosmology. 

For NFW halos, the value of the characteristic radius 
R s will depend on the cosmology used, the largest value 
obtained with the O^f = 0.3, — cosmology, fol- 
lowed by Om = 0.3, Qa — 0.7 and the Qm = 1, Oa = 
case. Now, it can be shown that for a fixed halo mass, 
M, the NFW potential will appro ach t he point-mass po- 
tential when lowering R s , see Fig. P15 Therefore, we ex- 
pect the closest resemblance to the point-mass case for 
= 1, Oa = 0, followed by Om = 0.3, Qa = 0.7 and 

n M = 0.3, n A = o. 

9. Limitations of method 

The method of HW, or rather the numerical implementa- 
tion of it, has some limitations due to the approximations 
used. Since the area of the light beam when exiting a cell 
is computed by linear interpolation from the value before 



entering the cell [see Eq. (23)], we need very small cells 
or very many cells to get robust results. Small cells means 
shorter interpolations while a large number of cells tends 
to average out the error in every interpolation. 




-0.4-0.2 0.2 0.4 
mag(filled beam) — 



0.6 0.8 1 
mag(lensed) 



In fact, in the limit R s 



0, $nfwM 



t(r). 



Fig. 8. Luminosity distributions for 10 000 sources at red- 
shift z = 1 in a Qjsi = 0.3, fl\ = universe. The magni- 
fication zero point is the filled-beam value. The full line 
corresponds to the point-mass case; the dashed line is the 
distribution for SIS halos, and the dotted line is the NFW 
case. 
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Fig. 9. Luminosity distributions for 10 000 sources at red- 
shift z = 1.5 in a J7m = 1> = universe. The magni- 
fication zero point is the fillcd-beam value. The full line 
corresponds to the point-mass case; the dashed line is the 
distribution for SIS halos, and the dotted line is the NFW 
case. 



Fig. 11. Luminosity distributions for 10 000 sources at 
rcdshift z = 1.5 in a f!ju = 0.3, Q\ = universe. The 
magnification zero point is the filled-bcam value. The full 
line corresponds to the point-mass case; the dashed line is 
the distribution for SIS halos, and the dotted line is the 
NFW case. 




-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 



rnag(f illed beam) — mag(lensed) 

Fig. 10. Luminosity distributions for 10 000 sources at 
rcdshift z = 1.5 in a Qm = 0.3, £!a = 0.7 universe. The 
magnification zero point is the fillcd-beam value. The full 
line corresponds to the point-mass case; the dashed line is 
the distribution for SIS halos, and the dotted line is the 
NFW case. 
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Fig. 12. Luminosity distributions for 10 000 sources at 
redshift z = 2 in a Om = 1, = universe. The mag- 
nification zero point is the filled-beam value. The full line 
corresponds to the point-mass case; the dashed line is the 
distribution for SIS halos, and the dotted line is the NFW 
case. 
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Fig. 13. Luminosity distributions for 10 000 sources at 
rcdshift z = 2 in a J7m = 0.3, Qa = 0.7 universe. The 
magnification zero point is the filled-beam value. The full 
line corresponds to the point-mass case; the dashed line is 
the distribution for SIS halos, and the dotted line is the 
NFW case. 
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Fig. 15. Lensing potentials normalized to M = 1. The 
dotted line corresponds to the potential of homogeneously 
distributed matter (filled-beam); the short-dashed line is 
a SIS halo with d = R c ; the long-dashed line is a SIS 
halo with d = 0.2R C , and the full line corresponds to the 
point-mass case. 
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rnag(f illed beam) — mag(lensed) 

Fig. 14. Luminosity distributions for 10 000 sources at 
rcdshift z = 2 in a Om = 0.3, CIa = universe. The 
magnification zero point is the fillcd-bcam value. The full 
line corresponds to the point-mass case; the dashed line is 
the distribution for SIS halos, and the dotted line is the 
NFW case. 
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Fig. 16. Lensing potentials normalized to M = 1. The 
dotted line corresponds to the potential of homogeneously 
distributed matter (filled-beam); the short-dashed line is 
a NFW halo with R s = R c \ the long-dashed line is a NFW 
halo with R s = 0.01i? c , and the full line corresponds to 
the point-mass case. 
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Fig. 17. Fractional mean error in area vs redshift. 



In order to quantify this effect we have made some tests 
with empty and filled beams using cell sizes obtained from 



the Schechter matter distribution. In Fig. 17, we have plot- 
ted the mean fractional error in area vs. redshift. To get 
a mean fractional error less than one percent, one would 
need to have z > 0.1. Note that this is the average error of 
individual measurements; the average value will be very 
close to the real value even at very small redshifts. For the 
cosmology used (£Im — 0.3, S1a = 0.7), a redshift of ~ 0.1 
corresponds to ~ 60 cells. 



10. Multiple imaging 

Another drawback with the HW method is that there is no 
way to keep track of multiple images from the same source. 
This poses no problem when the time delay between a pri- 
mary image and secondary images is big enough to allow 
observational separation of the images. However, if image 
separation is not observationally feasible, we have to com- 
pensate for multiple images. This will be the case when 
considering microlensing, see Sect. (10.3). 

A typical case of multiple imaging is shown in Fig. |l8|. 
Note that the secondary image has undergone a caustic 
crossing and that its parity has been reversed (this cor- 
responds to A < 0). Also, the secondary image is demag- 
nified compared to the primary image (remember that it 
is the relation between Sfl and the area that determines 
the magnification). Since the quantity 5il is more or less 
constant in our method, areas of secondary images often 
becomes very big, i.e., very dim. 




Fig. 18. Multiple imaging. 



For a point mass lens, the time delay between two im- 
ages with flux ratio r is given by 



At = 2i? Schw (l + z L )r(r), 
where 



T ( r ) = \(Vr- 



1 



lnr) 



(54) 



(55) 



For a lens at redshift z = 0.5 and a flux ratio of two, we 
get 



At ~ 2 ■ 10"" 



M 



(56) 



As an illustration, to get a time delay of the order of weeks, 
we need M ~ 10 n M Q . 

The image separation in the same case (assuming a 
source at z = 1) is given by 



A0 ~ 3 • 10" 



M 

Mo 



(57) 



i.e., a mass of M ~ 10 11 M Q will give an image separa- 
tion of ~ 1" . In the realistic case that a galaxy is causing 
the lensing, the surface brightness of the lens may make 
the detection of secondary images impossible (Porciani & 
Madau|i~998|). 



10.1. A point-mass universe 

Now, let us discuss what kind of information we can obtain 
from secondary images. In Fig. 19 we see what a source 
would look like in a universe filled with point-mass lenses. 
Since the deflection angle from a point mass can be very 
large, there will, apart from the primary image, be one 
secondary image next to every point mass lying in the 
vicinity of the line-of-sight to the source (not to mention 
all images with two or more caustic crossings)]^] Most of 
these images will be very dim, i.e. not observable. 

In our implementation of the HW method, we can hap- 
pen to follow any one of these. However, the vast majority 
of the flux will almost always be in the primary image, 
the only exception being the case where we have a point 
mass very close to the line-of-sight (see Fig. |o|). In this 
case, the primary and the brightest secondary image will 
have almost equal flux. If we want to consider the total 
flux from a source we need to correct for these images. 



The deflection angle is given by a = 2_Rgchw/f ■ 
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Of course, the correction would have to be larger in the 
case of Fig. pp| . Multiple-image corrections are discussed 
in Sects. 



10.3 and 10.4 below 



Note that if we happen to follow any of the weak sec- 
ondary images in Fig. we do not have any information 
on the flux of the primary image, i.e., we would not be able 
to discriminate between the case of Fig. |l9| and Fig. 
From this we conclude that weak secondary images does 
not give any valuable information concerning the total flux 
of a source. 
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Fig. 21. SIS universe. 
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Fig. 19. Point mass universe, first case. 
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Fig. 20. Point mass universe, second case. 



10.3. One-lens approximation for point masses 

Here, we will discuss reasonable corrections when images 
are non-separable. Since we have a very large number of 
images with two or more caustic crossings contributing to 
the flux of the total image, we can not correct for each one 
of them. Instead, we will make a first correction for the 
brightest secondary image using the approximation that 
lensing effects from a single lens (cell) is dominant. This 
is the so called one-lens approximation. Possible tests of 
the validity of this assumption is discussed in Sect. [n|. Its 
plausibility is also reinforced by Figs. 8 and 9 in HW. 

We can now use standard analytical expressions for 
one-lens systems for this correction: 



fi c 



|2 M -1| 



(58) 



In fact, this formula is valid also for secondary images 
but one then has to include the sign in /i, i.e., /i 2 < 0. 



However, Eq. (58) is only valid for the brightest secondary 



image, and one should therefore not apply it to very weak, 
secondary images. Also note, that to make the correction 
self-consistent, the magnification should be defined as the 
ratio of the expected area if the most prominent cell is 
empty and the real area. It is easy to see that the most 
important effect of this correction is to make the high 
magnification tail more pronounced. 



10.2. A SIS-galaxy universe 

Filling the universe with SIS-type galaxies, the situation 
would be quite different. No longer can we have very large 
deflection angles, and multiple images will therefore be 
less common.^ Since we are typically dealing with galaxy 
masses, images will be separable and we need only to con- 
sider secondary images bright enough to be observable. 
Since for (almost) every source, there will be at most one 
secondary image, it is a more straightforward procedure 
to associate to every secondary image a primary image. 

6 For SIS lenses, the deflection angle is given by a = 
in (v/c) 2 j£|/£, where v is the velocity dispersion of the lens and 
£ is the impact parameter, i.e., the magnitude of the deflection 
does not depend on £. 



10.4- One-lens approximation for SIS lenses 

It is as straightforward to compute the total magnification 
using the one-lens approximation with a SIS-lens. For a 
primary image we get 

/W = 2(/xi - 1). (59) 
For a secondary image we get 

/i corr = 2(1 - fi 2 ). (60) 

Since images are separable in the majority of cases, this 
is perhaps most interesting to use as a consistency check 
for the total luminosity. 

11. Consistency checks 

In order to validate our implementation of the method, 
we have performed a number of tests. 
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Fig. 22. Sum of primary areas divided by N to t for point 
masses (triangles) and NFW halos (crosses) in a VLm = 
0.3, Q\ = 0.7 universe. The number of primary images for 
each data point is Ni — 20 000. 

Filled and empty beam distances 

First - in order to get an independent check - we 
have compared results from our implementation of the 
HW method using empty cells and cells with a homo- 
geneous dust component with analytical empty-beam 
and fillcd-bca m res ults (using the Angsiz routine, see 
Kayser et al. 



1995) 



This has been done for a variety 
of cosmologies and the agreement is good even at very 
high z (less than 1 % discrepancy up to z = 10). 
Comparison with HW 

Second, we have reproduced most of the plots pre- 
sented by HW. We have also tested that primary im- 
ages satisfy j^— > 1, see Fig. |22| 
Average luminosities 

Imagine a source situated at the center of the sphere in 
the lower panel of Fig. || and observers sitting on the 
sphere at some constant redshift. Since the number of 
photons emitted by the source, as well as the area of 
the observers sphere is constant, i.e. independent of 
any details of the distribution of matter as long as the 
average matter density is the same, observers will on 
average receive the same flux, regardless of the mat- 
ter distribution. This average will be equal to the FL 
value, 



(61) 



When computing (i for a source, one has to consider 
the total flux if there are multiple images. If we con- 
sider primary images only, without any corrections, we 
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Fig. 23. Average magnification of primary images (nor- 
malized to filled-beam values) for point masses (triangles) 
and NFW halos (crosses) in an VLm = 0.3, J7a = 0.7 uni- 
verse. The number of primary images for each data point 
is N = 20 000. 

would expect to obtain a [i slightly lower than the FL 
value (since we lose some flux in secondary images) . 
Note that this does not mean that the average area 
should equal the FL area. Instead, 
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(62) 



i=i 



i.e., the FL area is the harmonic mean of the event 
areas, not the arithmetic (remember that (A) Harm < 
(-A)Arit)- Since we lose some flux in secondary images, 
we have (considering primary images only) 



A F L < (A) Harm < (A) Arit ■ 



(63) 



See Fig. |23| for a plot of the average magnification using 
point mass lenses and NFW lenses. The result for SIS 
halos is not plotted since the result is indistinguishable 
from the NFW result. Note that since the deviation 
from the filled-beam value is due to the absence of the 
flux contained in secondary images, Fig. |2^ indicates 
that the flux in secondary images is very low when 
using NFW or SIS halos, even at very high redshifts. 
The relationship between this test and the area test of 



HW is discussed in Appendix C 
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12. Analytical fitting formulas 



Since the numerical work involved in computing the re- 
sults presented in this paper is quite extensive, it may 
be convenient to use simple analytical fitting formulas for 
the probability distributions for the different cosmologies, 
which we now give. 

The probability distribution of the deviations from 
fillcd-beam magnitudes induced by SIS and NFW lenses, 
A =mag(filled beam)— mag(lensed), can be parametrized 
by an exponential function folded with a gaussian distri- 
bution: 



SIS 



P(A) oc 



-(A- 



dm' 



(64) 



Thus, for any combination of cosmological parame- 
ters, the corrections due to gravitational lensing can be 
described by three parameters: the mean of the gaussian 
distribution, /i, the standard deviation a, and the magnifi- 
cation tail, r. The expression in Eq. ( |64| ) can be rewritten 
so that (after normalization) a simpler integral is invoked: 



where 



r/<7-2(A-n) 



crfc 



r/a - (A - n) 



2t 



erfc(x) = -^= / e x dx' 

V7T 



(65) 



(66) 



In Fig. |24| the fits to the magnification distributions 
are shown for SIS lenses for source locations z = 1, 2 and 
3 and for three different cosmologies. In Fig. |2^ the same 
source locations and cosmologies are used but the lenses 
are of NFW type. 

The parametrization of Eqs. (Q) and (65) may be used 
to study the shape of the magnification distribution as a 
function of redshift and cosmological parameters, as shown 
in Figs. ^ and ^7|. 

The estimates of the cosmological parameters may be 
further refined and the understanding of the mass distri- 
bution in galaxy halos could be deduced by fitting the 
parameters r, a and fi to the residual magnitude distri- 
butions in future high-redshift supernova searches, such 
as the SNAPsat project or at NGST. The accuracy of 
this method will be discussed in a forthcoming paper 
(Bergstrom et al., in preparation). 

13. Discussion 



In this paper, the method of Holz & Wald ( |1998| ) has been 
generalized to include the effect of non- vanishing pressure. 
A motivation for this exercise is to use this method as part 
of a large package for a detailed simulation of most aspects 
high-redshift supernova determination of the geometry of 
the universe (Bergstrom et al., in preparation). The possi- 
bility to allow for various matter contributions is then de- 
sirable. Another improvement has been to allow for a more 
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Fig. 24. Fits of magnification distributions for SIS lenses 
using the analytical expression in Eq. 



realistic mass profile of the galaxy halos, the Navarro- 
Frenk- White distribution, where we have given analytical 
formulas for the geodesic deviation integrals. With this 
improvement, and by adding matter in our model uni- 
verses according to a Schechter function with parameters 
fixed by observation, we think we have achieved a suffi- 
cient degree of sophistication to make realistic predictions 
of the lensing effects of present and future deep supernova 
searches. 

Our results show, in agreement with previous work 



(Holz & Wald |1998[ Metcalf & Silk |1999| ) that the (un- 
likely) case of having most of the dark matter distributed 
in the form of point masses can be quite easily distin- 
guished from more realistic halo models, if a sufficient 
number of supernovae can be detected at large redshifts. 
It will be more difficult to distinguish between a singu- 
lar isothermal sphere model and the N-body results of 
Navarro, Frenk and White, unless f)jvf is large and one can 
measure supernovae out to z ~ 2 (see Fig. [l2]). However, 
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using the analytical expression in Eq. ([35]) ters in Eqs. (fii]) and (|6|) for SIS lenses. 



this also means that our method of determining the prop- 
erties of the lens population by matching the measured 
Schechter function with the required Qm gives quite ro- 
bust results for the predicted magnifications of Type la 
supernovae over a broad range of reasonable halo models. 
This will be helpful when making realistic predictions for 
future deep supernova searches. 
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Appendix A: : J a p for the Navarro- Frenk- White 
matter distribution 

In this Appendix, we give analytic expressions for the in- 
tegrals ( p5| ) needed for the evolution of the optical matrix 
for a bundle of light which passes near a mass distribu- 
tion of the NFW form. We assume a spherically symmetric 



matter distribution, the impact parameter is b, the charac- 
teristic radius of the NFW model is R s and the cell radius 
is R c . When R s ^ b, it follows that 



(J NFW ) X X = 2C 



b 2 



^R 2 c -b 2 {R 2 -2R 2 c -b 2 ) 
R 2 (R S + R C )(R 2 - b 2 ) 



In 




VRj^(2R 2 c +b 2 ) ( 

b 2 Rl { + R> 



2R s {R 2 -2b 2 ) 

' b 2 \R 2 -b 2 \y 2 



f(Jb,R a ,Rc 



(J NFW ) Y Y 



-2C 



R 2 (R S + R C ) 



(A.l) 
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/ iNFWa 



f(b,R s ,R c ) 



6^ ln l 1 + ^ 

2R s sgn{R s -b) , \ 



(J NFW ) y A - = 
In 




m ln 



Re 



R c + b 



3/2 p 2 



R 2 + 3{R c + bY 
3R 2 b 2 



(J 



NFW\Y 
) Y 



-2C 



R 2 C {R C + b) 




^W^¥(2R 2 c +b 2 ) ( Re 
b 2 R 3 n \ b 



b 2 (R c + b) 
The constant C is given by 

C = PcritScR^ = 



M 



(A.3) 



(A.4) 



where x c = R c /R s . 

Appendix B: : Parameter values for mass 
distribution 

Here, we have compiled various parameter values, relevant 
for the mass distribution discussed in Sect. |?]. In what 
follows, absolute magnitudes are denoted by M, while 
masses are denoted by M. 



Lin et al. See Eq. (32) of (Lin et al. |1996| ). 

• a = -0.70 ±0.05 

• M* - 5 log h = -20.29 ± 0.02 

• n* = (1.9 ± 0.1) • 10" 2 h 3 Mpc" 3 

Lin et al. gives its range of validity as —23.0 < M. < 
— 17.5, which is about ±2.7 magnitudes around their 
characteristic magnitude —20.29. This translates into 
a luminosity range L <E L„10 " 4 ' 2 ' 7 , and a mass range 
M G Af»10°- 4 - 2 - 7/(1 ~ /3) , which, with f3 = 0.2, gives 
M/M, e [0.045,22.4]. 



Kirshner et al. (Kirshner et al. 1983; KI) 
-1.25 (KI19) 
-1.25 (KI22) 
= -21.70 (KI18) 
= -22.70 (KI21) 
1.0 • 10 w h- 2 L Q (see Binney 



Tremaine 



(KI23) 



• aj = - 
ap 

• M*,j ■ 

• L*_ = 

pJ87| ). 

• n* = 1.2 • ! Q- 2 fr 3 Mpc -3 
Peebles (Peebles |1993[ PE) 

• /3 = 0.2 (PE3.39) 

• a = -1.07±0.05 (PE5.130) 

• M.* = —19.53 ± 0.25 + 5 log h in the bj magnitude 
system where M Q = 5.48. (PE 5.139, 5.140) 

• L» = 1.0 ■ lO lo e ±o - 23 /i _2 L (PE 5.141) 

• n* = 1.0 ■ lQ- 2 e ±QA h 3 Mpc~ 3 (PE 5.142) 
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• v* = 220 km/s (PE3.35) 

• A = 0.25 (PE3.35) 
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In the tanneries! work in tiite paper, we nave u«ed 

• /3 = 0.2 

• a = -0.7 

• i'tnin = 0.5-L* 

• imax = 2.0L* 

• n* = 1.9 • 10~ 2 /i 3 Mpc" 3 

• i>* = 220 km/s 

• A = 0.25 

Appendix C: : Mean area vs. luminosity 

Holz and Wald (|199S| ) use the fact that the area of the 
boundary of the past of an observer should be very nearly 
equal to the area of the past light cone in the underly- 
ing FL model. From this they deduce that adding up all 
areas - counting beams with an even number of caus- 
tics as positive and beams with an odd number as neg- 
ative - one should very nearly obtain the FL result, i.e., 
(Aq) ~ Afl, where the zero-subscript denote unweighted 
quantities. Here we investigate whether this claim is con- 
sistent with the fact that the average magnification should 
be equal to the FL magnification. 

Imagine that, without weighting, we get an area dis- 
tribution uq(A), where we assume that lensing effects are 
weak enough to only give primary images. According to 
HW, we then have 



ph/9908047 



astro- 



(Ao) 



Jn (A)AdA 
J n AdA 



A 



FL- 



(C.l) 



In our simulations, the area distribution will be a convo- 
lution of no(A) and a probability distribution p(A), 



n(A) =n (A)p(A) = n (A) 



A 



A 



empty 



(C.2) 



Since we want the harmonic mean of areas to be equal to 
the FL value, we get 



1 

JA) 



r 



n(A)/AdA 



J A -^n (A)/A cmpty dA _ 1 
J A -^ n (A)A/A empty dA (A )' 



(C.3) 



That is, the two consistency checks at least show consis- 
tency in the case where we only have primary images. The 
inclusion of non-primary images is highly non-trivial since 
the concept of negative areas is difficult to interpret as 
luminosities. One possibility is to combine this approach 
with the one-lens approximation method for multiple im- 
age correcting. However, this is difficult to do analytically 
since one would then need the form of no(A). 
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